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We analyze the strong hexagonal warping of the Dirac cone of Bi 2 Te 3 by angle-resolved photoe¬ 
mission. Along DM, the dispersion deviates from a linear behavior meaning that the Dirac cone 
is warped outwards and not inwards. We show that this introduces an anisotropy in the lifetime 
broadening of the topological surface state which is larger along FK. The result is not consistent 
with nesting. Based on the theoretically predicted behavior of the ground-state spin texture of a 
strongly warped Dirac cone, we propose spin-dependent scattering processes as explanation for the 
anisotropic scattering rates. These results could help paving the way for optimizing future spin- 
tronic devices using topological insulators and controlling surface-scattering processes via external 
gate voltages. 

PACS numbers: 73.20.-r, 73.20.At, 79.60.-i, 79.60.Bm 


Topological insulators (TIs) are characterized by an 
insulating bulk energy gap and gapless spin-polarized 
Dirac-cone surface states with elec tron spins locked per¬ 
pendicular to their linear momenta.li^ This peculiar spin 
texture is believed to pla y a central role in inducing exotic 
quantum phenoment!^, novel magnetic-spin physics,®^ 
as well as in the development of future spin-based low- 
power transistors among a variety of applications. Such 
perpendicular locking can be realized if spins of elec¬ 
trons occupying the Dirac cone move around in ide¬ 
ally circular constant-energy contours between the Dirac 
point and the Fermi level.lSl Owing to the spin-momentum 
locking, electrons on the surfaces of TIs are protected 
from backscattering,!^ an effect that might be of crucial 
importance in the generation of spin currents with re¬ 
duced dissipation,!^ coherent spin rotation,!^ spin-orbit 
qubits,!i^ as well as in other applications such as ma¬ 
nipulation of photon-polarization driven spin currents in 
real devices.^S! The spin texture of the Dirac cone can 
be affected by the presence of hexagonal warping,i.e., 
the Dirac cone is deformed in such a way that from the 
Dirac point to the energy of the Fermi level the constant- 
energy contours develop from an ideally circular shape 
to a hexagon, and subsequently to a snowflake-like Fermi 
surface. This strong distortion of the Dirac cone might 
open up new possibilities for observing other interesting 
phenomena, such as the breaking of time-reversal sym¬ 
metry and the surface quantum Hall effect under applied 
magnetic fields parallel to the surface,enhanced sur¬ 
face s catter ing or the formation of surface spin-density 

waves.fi^E^ 

In the context of exploring these possibilities, the 
observation and importance of hexagonal warping in 
the band dispersions of the topological surface states 


(TSSs),!i^lti^ as well as its theoretical description, 
has stimulated a manifold of theoretical and experimen¬ 
tal investigations on TIs. The strong influence of warping 
on the spin texture of TSSs giving rise to a finite out- 
of-plane spin orientation was theoretically predicted,!^ 
and by means of spin-resolved angle-resolved photoemis¬ 
sion (SR-ARPES) experimentally observed.!^ That in 
the presence of warping the electron spins are not any¬ 
more locked perpendicularly to their linear momentum 
was also proposed,!^ and the lifting of such perpendic¬ 
ular spin-momentum locking systematically investigated 
by means of SR-ARPES experim ents on TSSs of the pro¬ 
totypical TIs Bi2Te3 and Bi2Se3.!^^i2H2ll q^j^ggg investiga¬ 
tions revealed that in Bi2Te3, where warping effects are 
much stronger than in Bi2Se3,!^the distorted Dirac cone 
is characterized by a left-handed three-dimensional spin 
texture which features an out-of-plane spin component 
oscillating around the Fermi surface, and an in-plane spin 
component tangentially following the snowflake contour 
which results from strong warping.!^ 

In a similar way, the change that warping induces in 
the ARPES signal obtained from circular dichroism in 
the angular distribution was attributed to the lifting of 
the perpendicular spin-momentum locking,!^ or alterna¬ 
tively to the strong influence that warpi^ induces on 
the orbital angular momentum of TSSs.!^ The influ¬ 
ence of warping on the spin orientation was additionally 
predicted to alter the possible channels for quasiparti¬ 
cle scattering, and this effect has been the central topic 
of various theoretical investigations on TSSs taking into 
account spin-orbit scattering!^ or scattering at magnetic 
point defects.!^ The effect of quasiparticle interference in 
TIs was additionally studied by means of scanning tun¬ 
neling microscopy (STM) and spectroscopy (STS) and 
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from this type of experiments, an enhanced magnitude 
of the surface patterns resulting from electronic interfer¬ 
ence was observed and ascribed to hexagonal warping. 

Moreover, the possibility of a spin-density wave 
and Friedel-like oscillations of the local density of states 
was predicted^ and the latter experimentally observed 
by means of STM experiments.^ Equally important, a 
current-induced spin polarization was theoretically pre¬ 
dicted and it was shown that such a current may enhance 
the spin component perpendicular to t he su rface when 
significant hexagonal warping is present. ^23111 

In view of these observations and predictions, partic¬ 
ularly taking into account that warping can i nfluence 
the spin orientation of the surface-state electrons^^^HlIIll 
and that the spin orientation might a lter t he scatter¬ 
ing properties of the surface electrons, 1^3^ it is criti¬ 
cally important to especially investigate the momentum- 
dependence of the lifetime broadening of TSSs in the 
presence of strong hexagonal warping. Namely, whether 
a large Fermi-surface distortion is a weak perturbation on 
the lifetime broadening of TSSs or there are preferential 
momentum-space directions where this distortion results 
into a substantial modiflcation of the scattering rates of 
the surface-state electrons. ARPES represents an ideal 
technique for such a purpose, as it is a well-established 
method to determ ine the lifetime broadening of elec¬ 
tronic states,not only in strongly correlated systems 
such as transition- metal oxides,!^ 4/ rare-earths,^high- 
Tc superconductor^^nmior 3d ferromagnets,l23Hilbut also 
in Tls.l^smil Moreover, ARPES has recently become the 
most powerful tool in systematically revealing the linear 
dispersion of TSSs in energy and momentum space. 

In the present work, we utilize the TSS of the prototype 
TI Bi 2 Te 3 to analyze the strong hexagonal warping of the 
Dirac cone by angle-resolved photoemission. We investi¬ 
gate the impact of warping on the scattering properties 
of the surface-state electrons in this system. We demon¬ 
strate that the large distortion introduced by warping 
leads to an anisotropy in the lifetime broadening of the 
TSS which is larger along the FK direction of the surface 
Brillouin zone (SBZ). We further identify the underlying 
mechanism that gives rise to anisotropic scattering rates 
of the surface-state electrons. This result could help un¬ 
derstanding how to control surface-scattering processes 
that are relevant for spin injection in future spintronic 
devices based on TIs. 

We perform ARPES measurements at a temperature 
T =30 K using p + s linearly-polarized undulator radia¬ 
tion in ultrahigh vacuum better than 1 • 10“^° mbar with 
a Scienta R8000 electron analyzer at the UE112-PGM2a 
beamline of BESSY II. Bi 2 Te 3 single crystals are grown 
by the Bridgman method and cleaved in situ, with the 
sample temperature kept at 30 K. We use 21 eV and 
55 eV photons in combination with the sample geome¬ 
try shown in Eig. 1(a), where emitted photoelectrons are 
collected along the electron wave vector kg at an angle 
(j) = 45° with respect to the incident photon beam, which 
lies along the FK or FM directions of the SBZ within the 



FIG. 1: (Color online), (a) Geometry of the ARPES exper¬ 
iment in the laboratory reference frame. The light incidence 
and electron emission planes are indicated by black (dark) 
and red (light) dashed lines, respectively, (b) LEED pattern 
of the BbTea (111) surface acquired with an electron beam 
of 25 eV. 


yz plane of the laboratory reference frame. ARPES spec¬ 
tra are taken by rotating the sample about the x axis or 
z axis, while keeping the rotation about the y axis fixed. 
We use an analyzer entrance slit size of 0.2 mm x 25 
mm in width and length, respectively. Note that in the 
laboratory reference frame, ky and ky wave vectors 
are fixed parallel and perpendicular to the analyzer en¬ 
trance slit, respectively. The long axis of the entrance 
slit is parallel to the electron emission plane [red (light) 
dashed lines in Fig. 1(a)]. The angular resolution of the 
ARPES experiment parallel and perpendicular to the slit 
is about 0.1°, and the energy resolution ^30 meV. The 
high quality of the achieved (111) surfaces is verified by 
the sharp features in the low-energy electron diffraction 
(LEED) patterns [see Fig. 1(b)] and in angle-resolved 
photoemission of the valence band [see Fig. 2]. 

Figure 2(a) displays high-resolution ARPES disper¬ 
sions of the TSS, bulk conduction band (BOB) and bulk 
valence band (BVB) states of Bi 2 Te 3 sampled with 21 
eV photons around the F point of the SBZ along both 
in-plane ky ,,, and ky j, momenta. The differences in the 
electron dispersions for the two high-symmetry directions 
FM and FK sampled along ky are emphasized in Fig. 
2(b). Figure 2(c) displays the corresponding constant- 
energy surfaces extracted at different binding energies 
from the full photoemission mapping in Fig. 2(a). Fig¬ 
ures 2(d), 2(e), and 2(f) display analogous results to those 
obtained in Figs. 2(a), 2(b) and 2(c), respectively, but 
using a photon energy of hiz=55 eV. Quantization ef¬ 
fects due to band bendingl^SMl are absent, indicating 
that our ultrahigh vacuum experimental chamber con¬ 
tains an extremely small amount of residual gas which 
is not sufficient to generate a pronounced band bending 
effect giving rise to quantum-well states. A gapless Dirac 
cone representing the TSS with a Dirac point located at 
a binding energy (BE) of ^360 meV is clearly observed. 
The BE position of the BCB crossing the Fermi level, 
which at 21 eV reaches its minimum at a BE of ~35 meV 
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FIG. 2: (Color online), (a) High-resolution ARPES dispersions of the topological surface state, bulk conduction band and bulk 
valence band states of Bi 2 Te 3 sampled using 21 eV photons around the P point of the surface Brillouin zone along both in-plane 
k||_j. and momenta, (b) Comparison between the ARPES dispersions obtained along PM and PK. Both high-symmetry 
directions were oriented parallel to ky 2 ,. (c) Constant-energy surfaces extracted at different binding energies from the full 
photoemission mapping in (a). (d),(e),(f) Analogous results to those obtained in (a),(b),(c), respectively, but using a photon 
energy of hi/=55 eV. 


[see Fig. 2(c)], indicates that the crystals are intrinsi¬ 
cally n-doped. The intensity of the Dirac point strongly 
changes with the photon energy due to photoemission 
matrix elements,^ which also lead to changes in the rel¬ 
ative intensities from the TSS and bulk bands at differ¬ 
ent photon energies [compare, e.g.. Figs. 2(a) and 2(d)]. 
Moreover, because in ARPES the photon energy selects 
the component of the electron wave vector perpendicu¬ 
lar to the surface kj^, the contribution from bulk states 
to the ARPES intensity at a particular photon energy is 
also affected by their characteristic kj^dispersion. As a 
result of these effects, the BCB crossing the Fermi level, 
which clearly contributes at 21 eV, is strongly suppressed 
and virtually disappears at 55 eV photon energy [Figs. 
2(d)-2(f)]. In contrast, the energy-momentum dispersion 
of the TSS, despite the additional changes in its intensity 
distribution, remains clearly the same when varying the 
photon energy from 21 eV to 55 eV [compare, e.g., Figs. 
2(b) and 2(e)]. 

The TSS is strongly distorted by hexagonal warping. 


As a result, the constant-energy contours shown in Figs. 
2(c) and 2(f) progressively evolve from an ideal circular 
shape between the Dirac point and a BE of ^190 meV, to 
a hexagon up to a BE of ^110 meV, to progressively more 
pronounced snowflake-like constant-energy contours in a 
BE range of about 100 meV below the Fermi level. On 
the other hand, the lower Dirac cone, which emerges once 
the TSS bands disperse back to lower BE once they cross 
at the Dirac point, follows the energy-momentum disper¬ 
sion of the BVB, which at 21 eV reaches its maximum at 
a BE of ~230 meV along the FM direction of the SBZ 
[see Fig. 2(b)]. As a result, an overlap between the lower 
Dirac cone and BVB can be observed at this photon en¬ 
ergy. This is also the case at 55 eV, where the intensity 
from the BVB, although weaker, clearly contributes to 
the measured ARPES intensity, in contrast to the BCB 
which is strongly suppressed [see Figs. 2(d)-2(f)]. At 
21 eV photon energy, where the BCB intensity is very 
large, a strong overlap between BCB states and the TSS 
can be clearly distinguished along the FM direction in a 







4 


narrow BE range of ~20 meV below the Fermi surface 
[see, e.g., Figs. 2(a) and 2(c)]. This overlap is due to 
the fact that the BCB is also affected by the hexagonal 
distortion which follows the three-fold symmetry of the 
crystal structure around the [lllj-direction. In conse¬ 
quence, the BCB extends towards the TSS snowflake-like 
Fermi surface and touches it only along the FM directions 
of the SBZ. We note that the BVB also exhibits a dis¬ 
torted three-fold symmetric pattern but with somewhat 
different intensity distribution. In spite of this distor¬ 
tion, only the lower Dirac cone overlaps with the BVB, 
and this overlap seems to extend over a wide BE range, 
from around the Dirac point towards higher BE values. 

In order to quantitatively determine the changes in the 
lifetime broadening of the TSS in the presence of hexago¬ 
nal warping, in the following we perform a detailed anal¬ 
ysis of the TSS spectral width using momentum distri¬ 
bution curves (MDCs) extracted from the ARPES dis¬ 
persions shown in Figs. 2(b) and 2(e), where both TK 
and TM high-symmetry directions of the SBZ are ori¬ 
ented along k|| 2 , wave vectors (i.e., parallel to the ana¬ 
lyzer entrance slit). We use the standard method where 
the MDCs are fitted with Lorentzian peaks convoluted 
with a G aussian function representing the momentum 
resolution. ISHmU The results of our analysis using this 
procedure are shown in Fig. 3, and are only valid un¬ 
der the assumption that TSSs are weakly-interacting 
states. Note that in our evaluation the scattering rates 
of the surface-state electrons T = 211" = uoAk(i?) 
are directly related to the half width at half maximum 
(HWHM) of the Lorentzian peaks (Ak(F)/2), where 
E" =ImE(k, F) is the imaginary part of the complex 
self-energy S(k, A) =ReS(k, A) -|- iImE(k, F), and vg is 
the bare group velocity. This assumption is valid un¬ 
der the approximation in which vg is the renormalized 
group velocity Vg, implying that the real part of the self¬ 
energy E' is small, in agreement with recent observations 
on TSSs where a weak renormalization due to electron- 
phonon coupling in the vicinity of the Fermi level and 
slowl y vary ing energy-dependent values of E" have been 
found.l2Sll3 We p oint out that due to such a BE depen¬ 
dence of E",HSIi21the surface-state peaks still show a slight 
asymmetry with increasing BE, but by fitting MDCs with 
convoluted Lorentzian functions we have nevertheless ob¬ 
tained accurate results. 

Figures [3 (a) and [^(b) show few-selected fits [black 
(dark) solid lines] to experimental MDCs [red (light) dot¬ 
ted lines] obtained by varying the BE along TM and TK 
directions, respectively, and using 55 eV photons. Similar 
results were obtained at 21 eV by fitting the additional 
contribution from the BCB with an extra Lorentzian 
peak introduced into the analysis procedure, in similar 
way as previous reports at 8 eV photon energy for aged- 
Bi 2 Se 3 surfaces .1^ Figure [^(c) shows the corresponding 
surface-state ky 3 ;(F) dispersions extracted from the fit¬ 
ted peak BE positions along the TM and TK directions, 
which are represented by red (light) and blue (dark) solid 
lines, respectively. Along TK, the dispersion differs only 


slightly from a linear behavior. We have obtained accu¬ 
rate fits from the Fermi level up to BE values of ^260 
meV. In the case of the lower Dirac cone, the above 
mentioned contribution from the BVB at 21 and 55 eV 
prevents an accurate fitting. The different behavior of 
the fitted dispersions along the two high-symmetry di¬ 
rections, as seen in Fig. K c), is in agreement with the 
presence of hexagonal warping. However, our experimen¬ 
tal results are not consistent with the expectations oi k-p 
theory, which predicts a strong deviation from the linear 
dispersion with increased group velocity along the TK di¬ 
rection.In reality, we observe a deviation to reduced 
group velocities along TM, meaning that the Dirac cone 
is warped outwards and not inwards along this direction. 
Such a disagreement between k -p theory and experiment 
is due to the theoretical energy-momentum dispersion, 
despite the good agreement with the experimental Fermi- 
energy contour. The origin of the discrepancy comes from 
the effect of bulk conduction bands, which come down in 
energy along the TM direction. The level repulsion be¬ 
tween these bulk bands and surface states causes surface 
bands bending down in TM direction, leading to the ob¬ 
served hexagonal Fermi surface. The surface-only k ■ p 
Hamiltonian cannot capture this interaction with bulk 
bands, in contrast to first-principle calculations.^ 

The anisotropic behavior of the surface-state disper¬ 
sions is also reflected in the HWHM of the Lorentzian 
fits, as shown in Fig. |^(d) for 21 and 55 eV photons, 
which provide similar results (solid and dotted lines, re¬ 
spectively). In the BE range where the dispersion along 
the TM direction starts to deviate from the linear be¬ 
havior, i.e., from about 210meV BE in the present case, 
the HWHM increases from a value of ~0.01lA“^ to a 
maximum value of ~0.015 at a BE of ~80meV and 
then decreases towards the Fermi energy to approach the 
value of ~0.009A“^, which approximately matches the 
one along the TK direction. If we would assume a con¬ 
stant and isotropic bare group velocity vg for both di¬ 
rections, this would mean enhanced scattering for elec¬ 
trons moving along TM as compared to electrons moving 
along TK. However, the experimental group velocities 
are anisotropic and not constant as shown Fig. [|(e). We 
do note that using the definition of the group velocity 
as dE/hd\i, implies that it has to be determined from 
the A(k) dispersion, and not by differentiating the k(A) 
dispersion and then calculating the reciprocal value of 
the result. This method would be only valid for a lin¬ 
ear dispersion where the group velocity is nearly BE- 
independent. In the case of Bi 2 Te 3 , one has either to ap¬ 
proximate the group velocities linearly over small energy 
intervals, taking into account artificial energy steps, or to 
determine the £'(k) dispersion. For the results shown in 
Fig. |( e), we have chosen the latter method and mapped 
the function Ug(k) into Vg{E). To avoid the introduction 
of additional noise, we have fitted the resulting functions 
with polynomials. As seen in Fig. [^(e), along the TK 
direction the group velocity varies only slowly with BE 
and, in particular, it increases towards lower BE follow- 
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FIG. 3: (color online). Analysis of MDGs of the surface state of Bi 2 Te 3 . (a),(b) Few-selected fits [black (dark) solid lines] to 
experimental MDGs [red (light) dotted lines] obtained by varying the BE along (a) FM and (b) FK directions, respectively, 
and using 55 eV photons, (c) The k||^ 2 .(i?) dispersion of the TSS along the FM and FK directions [red (light) and blue (dark) 
solid lines, respectively] shows an anisotropic behavior, (d) The anisotropy is also reflected in the HWHM of the Lorentzian 
peaks, which for 21 and 55 eV photons provide similar results (solid and dotted lines, respectively), (e) The group velocity, 
obtained from a polynomial fit (solid lines) to the derivative of the i?(k) dispersion (open circles), is strongly altered over the 
analyzed BE range along the FM direction (eVA= hm/s). (f) Multiplication of the HWHM with the group velocity at each 
BE gives the imaginary part of the self-energy E", which reflects scattering rates that are anisotropic in k-space. 


ing a linear dependence, in agreement with the slight 
parabolic shape of the TSS dispersion between the Dirac 
point and the Fermi level. In contrast, along the FM 
direction, we find a more pronounced BE dependence of 
the group velocity which varies between ~3.1eVA and 
the minimum value of ~1.5eVA reached around a BE of 
^80 meV. Specifically, the Eermi velocity Ui? is different 
along the two high-symmetry directions, being ^2.6 eVA 
and ^3.5 eVA along FM and FK, respectively. 

If we now use the standard method of multiplying the 
HWHM values extracted from the fits to MDGs with the 
experimental group velocities to determine the imaginary 
part of the self-energy E", we obtain the results shown in 
Fig.|(f). We do note that these results are qualitatively 
reversed as compared to the ones shown in Fig.[^(d). Af¬ 
ter multiplication by the group velocities, we find scat¬ 
tering rates which are almost constant in BE along FK, 
with a slight and smooth decrease while approaching the 
Fermi level. At a BE higher than ^200 meV, the scatter¬ 
ing rates are qualitatively similar for both FK and FM 
directions, as expected for a circular Dirac cone. How¬ 


ever, the scattering rates for electrons moving along FM 
are strongly reduced as the dispersion experiences more 
and more the hexagonal deformation. From a BE of ~200 
meV towards the Eermi level, an almost linear decrease 
of the scattering rates is clearly observed along this di¬ 
rection. Below a BE of ~100 meV, the linear decrease 
slowly saturates, and nearly constant values of E" are 
reached down to a BE of ~25 meV. In the narrow BE win¬ 
dow between ^25meV and the Fermi level, E" slightly 
increases to approach the scattering rates along FK, ex¬ 
hibiting a significant upturn behavior in the immediate 
vicinity of the Fermi level. We point out that in this BE 
region, a pronounced overlap between the TSS and BOB 
is clearly observed along FM with 21 eV photons, as dis¬ 
cussed above [see Figs. 2(a) and 2(c)]. One would expect 
that at 21 eV photon energy, the intensity contribution 
from the BCB in the overlapping region prevents us from 
fitting more reliably the TSS peaks in this small BE en¬ 
ergy window. However, for 55 eV photons, the BCB is 
strongly suppressed, and the upturn behavior of E" along 
FM can be clearly observed in Fig.[^(f). This finding al- 
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FIG. 4: (Color online). Extraction of undistorted bands from the constant-energy surface mapping at 21 eV. (a) From cuts 
perpendicular to the high-symmetry directions a sequence of i5(k||_2;) (£'(k||_y)) dispersions at different constant k||_y (k||,a;) 
values is obtained (see text), (b) EDCs extracted from the sequence at k||^j,=0 (k||^j,=0) along the FM (FK) direction. The 
EDCs cut through a flat band and show rather symmetric peaks for the surface state, as represented by colored squares. Solid 
black (dark) lines show the result of the fits, (c) The imaginary part of the self-energy E” shows the same qualitative behavior 
as the results obtained from MDCs in Fig.[^(f). Solid and dashed lines are results for 21 and 55 eV photons, (d) The differences 
in the lifetime broadening of the TSS can also be noticed by directly comparing the width of the TSS peaks extracted at a BE 
of ~80 meV [marked by a vertical dashed-gray (light) line in Fig.[^(c)]. Results obtained with p -\- s linearly (filled symbols) 
and positive circularly polarized light (open symbols) are shown. In (c) and (d), blue (dark) and red (light) colors denote F-K 
and F-M directions, respectively. 


lows us to attribute the upturn behavior very near the 
Fermi level to an increased surface-bulk coupling within 
the overlapping region, which leads to an enhancement 
of the scattering rate of the surface-state electrons along 
FM. As a result, the scattering rates along FK and FM 
directions behave more isotropically at the Fermi level. 
It should be emphasized that this type of coupling is in¬ 
dependent of the photon energy, despite the fact that 
the spectral weight of the BOB is strongly suppressed for 
55 eV photons. In this respect, we would like to point 
out the basic fact that the surface-projected bulk band 
structure, which hosts the signatures of the coupling, is 
independent of the probing photon energy (effectively 
kj_ values). In other words, if the overlap between TSS 
and BCB along FM direction is observed at a particular 
photon energy (21 eV in our case), the signatures of en¬ 


hanced surface-bulk coupling in the linewidth of the TSS 
will remain photon-energy independent. This will be the 
case even if the BCB disperses with kj_or its intensity 
is suppressed at a particular excitation energy (such as 
55 eV in our case). Similarly important, it should be 
noted that around 80 meV BE, below the BCB mini¬ 
mum, the anisotropy between the scattering rates along 
FK and FM is rather large for both 21 and 55 eV photons. 
Around this BE, the corresponding constant-energy con¬ 
tours already exhibit a snowflake-like shape with a level 
of distortion similar to the one of the Fermi surface [see 
Figs. 2(c) and 2(f)]. This result strongly indicates that 
away from the immediate vicinity of the Fermi level, the 
anisotropic behavior of E" observed in Fig.|^(f) is related 
to warping. 

We emphasize that the results shown in Fig.[^(f) have 
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not been corrected for impurity-scattering contributions, 
which in a trivial way increase E" by a small constant 
and energy-independent value represented by the offset 
at the Fermi level. Other extra contributions, such as 
the one from electron-phonon coupling, can be found in 
the Debye model for the high temperature limit Here 
the electron-phonon broadening depends linearly on the 
temperature T as E"_p^=7rAfcB'F, where ks and A are 
the Boltzmann and electron-phonon coupling constants. 
From temperature-dependent measurements in our own 
Bi 2 Te 3 samples we find a weak electron-phonon coupling 
characterized by A « 0.17, in agreement with recent 
studies.l^This result leads to « 2 meV, which rep¬ 

resents an almost negligible energy-independent contri¬ 
bution to the lifetime broadening due to electron-phonon 
interaction. Other contributions to the linewidths due 
to final-state broadening, which for two-dimensional elec¬ 
tronic states are re duced to an extra broadening of purely 
geometrical origin, are also negligible. In this case 
the measured peak widths are compressed or expanded 
depending on a geometrical factor C so that the mea¬ 
sured lifetime broadening is represented by CS", being 
C = 1/11 — TOUosin^0|, Vo the group velocity, m the elec¬ 
tron mass and 9 the corresponding emission angle.^^^^^ 
Due to the small size of the Fermi surface of the TSS, 
our detail analysis reveals that for 21 and 55 eV photons 
C « 1 along both FM and FK directions. Therefore, 
our results of the MDCs analysis show the importance 
of a proper group velocity determination and mainly de¬ 
pend on its specific behavior, in particular taking into 
account the obvious differences between the HWHM and 
the S" values obtained in Figs.|^(d) and [^(f). Hence, 
in order to confirm the validity of these results and, be¬ 
fore discussing the behavior of E" in more detail, in the 
following we want to justify the usage of the experimen¬ 
tal group velocities by a different method, the results of 
which are shown in Fig.|^ 

Differently from the method used to obtain the results 
of Fig.[^ in Fig.|^ we perform an analysis of the TSS 
spectral width using energy distribution curves (EDCs) 
extracted from the ARPES dispersions shown in Eigs. 
2(a) and 2(d), where the FK and FM directions are ori¬ 
ented along k|| 2 ; and ky ^ wave vectors, respectively (i.e., 
parallel and perpenducular to the analyzer entrance slit). 
This method allows us, on the one hand, to avoid an ar¬ 
tificial broadening caused by the angle between the band 
dispersion and the cut direction and, on the other hand, 
to check whether the sample geometry influences the re¬ 
sulting scattering rates by directly comparing to the re¬ 
sults of Fig.[^ Moreover, if the self-energy E(k,E) is 
slowly varying with BE, the HWHM of the TSS peaks in 
each EDC gives the imaginary part of the self-energy E". 
This method has been used, for instance, to study many- 
body effects in Mo(llO) surfaces,^3d ferromagnets,l32lEil 
or high-Tc superconductors.!^ 

The procedure we have used is illustrated in detail in 
Eig.[^ At the top-left side of Fig.[^(a), as a reference, 
we show a constant-energy surface extracted at a BE of 


~80meV from the 21 eV data of Eig. 2. In order to 
quantitatively analyze the peak widths corresponding to 
the TSS bands along both FM and FK directions using 
EDCs, we extract two sequences of A(kj[) plots from the 
21 eV volumetric data set shown in Fig. 2(a). These 
sequences cut perpendicularly to the FM and FK direc¬ 
tions, as indicated on top of the constant-energy surface 
of Fig.|^(a) by dashed-red (light) and blue (dark) lines 
oriented parallel to ky and ky ^ wave vectors, respec¬ 
tively. Few-selected results of the A(ky) plots belonging 
to each sequence are shown on the right and bottom sides 
of Fig. ^(a). Note that along FM, the plots are displayed 
from the top to the bottom of Fig.[^(a), while along FK, 
from left to right, which corresponds to increasing values 
of ky y and ky j,, respectively. In each of these plots, the 
TSS band reaches a minimum which is represented by a 
small energy region located at ky_3;=0 (ky_y=0) along FM 
(FK). Clearly, these band minima change in BE position 
through the sequences. In order to obtain the imaginary 
part of the self-energy E", we analyze the peak widths 
of these band minima along FK (FM) using the corre¬ 
sponding EDCs extracted at ky ^=0 (ky 3,=0) [marked 
by dashed-white (light) lines in each E(k|[) plot]. This 
way, we ensure that the broadening of the peaks does 
not suffer from the aforementioned artificial broaden¬ 
ing. Moreover, because the geometrical factor C is « 1 
along both FM and FK directions, we do not expect ex¬ 
tra broadening of purely geometrical origin contributing 
to the EDC widths. In Fig.[^(b), we show few-selected 
fits to EDCs extracted for 21 eV photons from the se¬ 
quence at ky_2;=0 (ky_y=0) along the FM (FK) direction 
[top (bottom) panel]. The experimental EDCs, which are 
represented by colored-dotted lines, exhibit rather sym¬ 
metric peaks following the energy-momentum dispersion 
of the TSS. Each color corresponds to a different cut 
through the sequences. In analogy to the method used 
in Fig.|^ the peaks are fitted by Lorentzian profiles con¬ 
voluted with a Gaussian function representing the energy 
resolution (see, for instance, Ref. |33]for more details on 
the procedure). Similar analysis was performed for the 
55 eV data of Fig. 2(d). Although we find almost absent 
asymmetries caused by the BE dependence of E", other 
problems connected to the analysis of the EDCs might 
still exist, like the secondary electron background or de¬ 
viations near the cut-off at the Fermi level. Nevertheless, 
in our case we find a low-background signal and rather 
symmetric peaks across the Fermi level, in a similar fash¬ 
ion as in previous studies. 

The E" values, which are directly represented by the 
HWHM of the Lorentzian peaks obtained by fitting the 
corresponding EDCs, are shown in Fig.[^(c) for both FK 
[blue (dark) circles] and FM [red (light) circles] high- 
symmetry directions. The analysis for 21 and 55 eV pho¬ 
tons provides similar results (solid and dotted lines, re¬ 
spectively). Again, pronounced differences in the behav¬ 
ior of E" along FM and FK leading to a remarkable mo¬ 
mentum dependence of the scattering rates, are clearly 
observed. Such differences in the lifetime broadening of 







the TSS are independent of the light polarization, as it 
can also be noticed in Fig.|^(d) where we directly com¬ 
pare for linearly and circularly polarized 21 eV photons of 
positice helicity the width of the TSS peaks extracted at a 
BE of ^80 meV [marked by a vertical dashed-gray (light) 
line in Fig. El( c)]. Note that this BE corresponds to the 
constant-energy contour shown in the top-left side of 
Fig.ga), where the hexagonal distortion of the constant- 
energy surfaces and the anisotropy of the scattering rates 
are rather large. We point out that similar results are ob¬ 
tained either by reversing the light helicity of the circular 
polarization to negative or using s-polarized light. Our 
results are consistent with unobservable differences be¬ 
tween the orbital-projected scattering rates and with the 
fact that the TSS signal within the first few atomic lay¬ 
ers is dominated by a strong contribution from Pz-type 
orbitals. Moreover, comparing Eigs. if) and |^(c), we 
find qualitative agreement, a fact that corroborates our 
initial observation of reduced scattering rates along FM 
as compared to FK direction in the presence of hexago¬ 
nal warping. Hence, it is reasonable to assume that for 
the TSS of Bi 2 Te 3 , multiplication of the HWHM values 
obtained from MDCs with the experimental group veloc¬ 
ities is an equally valid method as the EDC analysis, at 
least to qualitatively look at the behavior of E" over a 
certain BE range. In this respect, we have to stress that 
the validity of both methods hinges on the assumption 
that the TSSs are weakly-interacting states. This fact is 
also evidenced in the EDC analysis by the excellent ac¬ 
curacy of the EDC fits using symmetric functions or, in 
other words, by the absence of asymmetries in the TSS 
peaks due to damping of quasiparticle excitations. Addi¬ 
tionally, the qualitative agreement between Eigs.|^f) and 
El( c) implies that our results are independent of the sam¬ 
ple geometry, indicating that other effects such as orbital 
spin interference^ do not play a role in our observation 
of anisotropic scattering rates. Despite the fact that very 
near the Fermi level the EDC analysis becomes more diffi¬ 
cult due to the peak intensity cut-off by the Eermi energy, 
it is also important to observe that along FM direction, 
the results of Fig. c) reveal a significant upturn behav¬ 
ior of E" in the immediate vicinity of the Fermi level. 
This observation is in qualitative agreement to the re¬ 
sults of the MDC analysis shown in Fig.|^(f), indicating 
that the upturn behavior and its independence on the 
photon energy are not connected to the sample geome¬ 
try or to the analysis method, but rather to surface-bulk 
coupling, as discussed above. 

In the following, we would like to focus on the main 
outcome of the present work which is the anisotropic 
effect of warping on the scattering rates in the region 
where the TSS and the BCB are decoupled. In par¬ 
ticular, we would like to examine various mechanisms 
that can explain this observation. These mechanisms, 
which are all related to warping and in some degree con¬ 
nected to each other, contribute with different weights, 
as we will discuss below. An important mechanism con¬ 
tributing to the anisotropic behavior of the scattering 


rates is spin-dependent scattering, i.e., spin-dependent 
decay of photoholes in our case. For a two-dimensional, 
non-degenerated, time-reversal-symmetric surface state, 
backscattering is forbidden, and the presence of strong 
warping does not affect this prohibition. However, this 
scenario is only valid for 180° backscattering, as scatter¬ 
ing under different angles is not completely suppressed. 
In particular, the probability for a scattering event away 
from 180° will be enhanced if there is a strong overlap 
between parallel spin configurations in the initial and fi¬ 
nal states involved in every scattering event.^ Equally 
important, as we have already mentioned, is the theo¬ 
retical prediction that warping affects the spin such that 
it is not anymore locked perpendicular to the electron 
momentum, but tangential to the constant-energy con¬ 
tours.Such predicted ground-state electron spin tex¬ 
ture is illustrated in Eig. 5 for a snowflake-like [Eig.l^a)] , 
a hexagon-shaped [Fig. 5(b)], and a circular [Fig^(c)[ 
constant-energy contour. Note that due to the limited 
momentum resolution of presently available SR-ARPES 
setups, a current experimental challenge is the acces¬ 
sibility to the pronounced tangential in-plane spin ro¬ 
tations seen in Fig.[^(a) near the FM direction. This 
type of tangential spin rotations have been claimed to be 
found experimentally in Bi 2 Se 3 ,l^ where the Dirac cone 
is only slightly distorted by warping, but the method 
used in these experiments, i.e., circular dichroism, ap¬ 
pears ques tiona ble as it reverses several times with pho¬ 
ton energy.l22llsl 

In the context of SR-ARPES, it should be emphasized 
that the recently observed changes in the spin polar¬ 
ization of photoelectrons emitted from TSSs as a func¬ 
tion of light polarization are due to the interplay be¬ 
tween spin-dependent matrix-elements, dipole-s election 
rules and photoelectron spin-interference effects.^^^EHIIl 
In other words, the observed effect of photoelectron 
spin manipulation with light polarization is not caused 
by ligh t-induce d changes in the TSS initial-state spin 
texture.l^^l^^^^ Moreover, our observation in Fig.|^(d) 
that the anisotropy of the scattering rates does not 
change with the photon polarization demonstrates that 
it is completely independent of light-induced changes in 
the photoelectron spin. Therefore, in the following we 
would like to focus our attention in the theoretically pre¬ 
dicted gro und-st ate spin texture of a strongly warped 
Dirac cone.*I^*^ This peculiar spin texture allows us to 
explain in detail our observation of anisotropic scattering 
rates in Bi 2 Te 3 via spin-dependent scattering, as we will 
discuss below. 

Let us assume that an electron has been removed from 
the vicinity of a tip of the snowflake-shaped constant- 
energy surface [Eig.l^(a)] along the FM direction. In this 
case, the electron spin will be fully parallel to the surface 
plane (i.e., Sx ^ 0 and/or Sy ^ 0, while = 0), consis¬ 
tent with the mirror symmetry of the crystal.^ Scatter¬ 
ing under small angles is allowed, but it is less likely to 
occur along the FM direction [e.g., qi in Fig.[^(a)], due 
to the fact that the spin varies strongly while following 
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FIG. 5: (Color online), (a)-(c) Schematics of different scattering channels, (a) For a snowflake-like constant-energy contour, 
the in-plane spin directions [components Sx and Sy, as indicated by blue (light) arrows] change stronger in the vicinity of the 
FM direction as compared to the FK direction, leading to a reduced probability of small-angle scattering events, e.g., qi as 
compared to events exemplified by q 2 . The out-of-plane spin rotation for k vectors away from the FM direction is sixfold 
symmetric, as indicated by the symbols © (sz points out of the paper plane) and © (sz points into the paper plane). The 
scattering probability is increased for scattering under qs, whereas backscattering (e.g., ^ 4 ) is still forbidden. The almost flat 
regions of the constant-energy contours, connected by red (light) scattering vectors {qneat) fulfil a near nesting condition, (b) 
Even stronger nesting is expected for a hexagonal shape, however, for this case time-reversal symmetry strongly suppresses 
scattering. Thus a spin-density wave would be favored over a charge-density wave (see text), (c) Similar scheme for a circular 
shape where scattering vectors such as 52 and ga contribute isotropically to the total scattering rate, (d)-(f) Simplified model 
calculations of the spin-dependent scattering amplitude, which includes contributions from all possible q vectors within (d) a 
snowflake-like, (e) a hexagon-shaped and (f) a circular constant-energy contour. 


the constant-energy con tour [as indicated by blue (dark) 
arrows in Fig.j^. In contrast, a hole caused by removing 
an electron from the same BE with its linear momentum 
pointing along the FK direction can be filled with elec¬ 
trons from adjacent sites in k-space, since the in-plane 
spin components are nearly parallel to each other [e.g., 
52 in Fig.[^(a)]. Furthermore, along the FK direction, 
the electron spin exhibits a finite out-of-plane compo¬ 
nent Sz [denoted by © and © symbols in Figs. IDla) and 
[^b)]. This spin component, which can lead to out-of- 
plane spin p olarization values expected to be as large as 
reverses its sign with an angular periodicity 
of ^ and vanishes along the FM direction due to mirror 
symmetry.Note that Sz varies only smoothly around 
the distorted constant-energy contours of Figs.j^a) and 
ib), meaning that even though adjacent sites in k-space 
might have different absolute values of the out-of-plane 
spin polarization, these values will be always of the same 
sign. This specific behavior of Sz strongly favours spin- 


dependent scattering as explanation for the increased 
scattering rates along FK direction. Such an argument 
can be further exemplified by considering a scattering 
vector like 53 in Fig.[^(a). The initial state (marked by 
the black (dark) dot at the beginning of 53 ) has finite spin 
components Sx and Sy, but there is no overlap between 
these spin components and the ones in the final state, 
which has Sy=0 and Sx pointing antiparallel with respect 
to Sx in the initial state. However, Sz is pointing out of 
the surface plane for both initial and final states, lead¬ 
ing to a relatively strong contribution to the scattering 
probability due to the large overlap between these two 
parallel out-of-plane spin configurations.^ In contrast, a 
comparable scattering event for electrons moving along 
the FM direction, such as the one exemplified in Fig.[^(a) 
by 55 , where Sz= 0 , will have only a small overlap be¬ 
tween the corresponding Sy in-plane spin components in 
both initial and final states. Similarly, in a hexagon¬ 
shaped constant-energy contour there will be always a 
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small overlap between parallel spin components for scat¬ 
tering vectors such as < 71 , (72 or <73 [see Fig. [^b)]. In the 
case of a circular constant-energy contour, because 8 ^= 0 , 
only contributions from scattering vectors such as (72 or 
(73 would be possible [see Fig. [11(c)]. 

In order to gain more insight on the overlap between 
different spin configurations, we have evaluated on the 
basis of a simplified model spin-dependent scattering ma¬ 
trix elements of the form | (s(k) | s(k-|-q)) |^. These matrix 
elements represent the scattering amplitude between two 
different spin states s(k) and s(k-|-q). The contribution 
to the total scattering is represented in our calculation 
by the sum of individual matrix elements for all possible 
q vectors around the constant-energy contour, weighted 
to the total spin polarization. Here s denotes the elec¬ 
tron initial-state spin configuration in three dimensions, 
which is represented by the predicted gr ound- state spin 
texture of a strongly warped Dirac cone.li^*^ This type 
of spin-selective scattering processes are enhanced for 
aligned spins and completely suppressed for opposite and 
perpendicular spin orientations.^ In our simplified model 
the constant-energy contours are modeled by paramet¬ 
ric equations, and the scattering amplitude in each k- 
space site is quantified by the inner product of three- 
dimensional spin vectors. Scattering events with zero 
or anti-parallel spin projections are strictly forbidden. 
The spin texture is represented by spin-polarization vec¬ 
tors with in-plane spin components tangential to the 
constant-energy contours. The magnitude of the out-of 
plane spin polarization along the FK dire ction is hxed 
to the theoretically expected value of and in¬ 

terpolated all around the constant-energy contour using 
a sinusoidal function with a periodicity of 27r/3. Simi¬ 
larly, we consider a constant magnitude of 50% for the 
expected value of th e total initial-state spin polariza¬ 
tion of the TSS.li^l^ Other effects, such as orbital spin 
interference,!^ are not taken into account. 

The results of the model calculations are depicted 
in Figs.[^(d), [^(e), and [^(f) for a snowflake-like, a 
hexagon-shaped and a circular constant-energy contour, 
respectively. Our simulation reveals that in the case of 
an ideal circle, spin-dependent scattering events are likely 
to occur isotropically. When the influence of warping de¬ 
forms the constant-energy contour from a hexagon to a 
snowflake-like shape, the anisotropy in the scattering is 
manifested as higher and lower spin-dependent scattering 
amplitudes along the FK and FM directions, respectively. 
These results are in qualitative agreement with our ex¬ 
perimental finding of anisotropic scattering rates, indi¬ 
cating that spin-selective scattering is the driving force 
behind this observation. We point out that including 
in our model unpolarized BCB states which only over¬ 
lap with the TSS bands along the FM direction, leads to 
a nearly isotropic scattering behavior. This is the case 
if we assume that enhanced surface-bulk coupling even¬ 
tually quenches the spin polarization in the overlapping 
region. Such a strong reduction of the spin polarization 
along FM opens up multiple scattering channels along all 


directions. This limiting case of our simplihed model is 
in principle consistent with the upturn behavior of E" 
along the FM direction. 

Another mechanism that might play an important role 
is nesting.In contrast to an ideally circular Dirac cone, 
or for example a two-dimensional free-electron gas, the 
warped constant-energy contours of Bi 2 Te 3 in the vicin¬ 
ity of the FK direction fulfil a nesting condition, espe¬ 
cially those of hexagonal shape, such as the one shown in 
Fig.[5](b). Note that for a snowflake-like constant-energy 
contour, at least a near-nesting condition must also be 
fulfilled.Nesting means that parts of the constant- 
energy contour can be mapped onto each other by a sin¬ 
gle scattering vector quest [indicated by the red (light) 
arrows connecting opposite sides of the constant-energy 
contours in Fig.[^(a) and [^(b)]. It is expected that ful¬ 
filling the nesting leads to enhanced scattering, but also 
to instabilities, like, for example, a charge-density wave. 
However, time-reversal symmetry disfavours the forma¬ 
tion of a charge-density wave,E 21 because independently 
of the warping, states at opposite ky momenta have op¬ 
posite spin [see, e.g., the blue (dark) arrows in Fig.[^(b)]. 
Alternatively, the formation of a spin-density wave in 
Bi 2 Te 3 has been also suggest ed as a possible instabil¬ 
ity introduced by the warping. Although we do not 
observe clear signs of a spin-density wave in the ARPES 
band dispersions of Fig. 2, i.e., a band gap at the Fermi 
level, this fact does not completely exclude a contribu¬ 
tion from nesting in the enhancement of the scattering 
rates along the FK direction. However, since the nesting 
vector quest given in Fig.[^(b) is strongly suppressed by 
time-reversal symmetry and the one in FigJ^(a) does not 
favour scattering of electrons along the FK direction, it 
is likely that nesting plays minor role in our observation 
of anisotropic scattering rates. 

Finally, scattering of photoholes through transitions 
into the bulk continuum might additionally contribute 
to the scattering rates because of the presence of dis¬ 
persing bulk states which do not directly overlap with 
the TSS bands. In an analysis of the scattering rates of 
the surface-state electrons of Bi 2 Se 3 ,!^ a pronounced de¬ 
crease of S" was found from the minimum of the BCB 
towards the Fermi level along the FK direction. This ob¬ 
servation was attributed to scattering processes strongly 
dominated by transitions into bulk electronic states, a 
finding that is only in qualitative agreement with the 
overall behaviour of our data for Bi 2 Te 3 along the FK di¬ 
rection. Note that we observe average values of E" that 
change very smoothly and by a small amount between 
the minimum of the BCB and the Fermi level along this 
direction. Moreover, the upturn behaviour of E" in the 
immediate vicinity of the Fermi level along the FM direc¬ 
tion shows the opposite trend that the one expected from 
this type of scattering channels. Likewise, in a BE range 
around the Dirac point, parts of the TSS are surrounded 
by the BVB and photoholes can also be scattered through 
transitions into the bulk continuum. As a result, a strong 
decrease of the scattering rates would be equally expected 
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from high BE up to around the BVB maximum. How¬ 
ever, our results reveal that in this BE region the changes 
in E" are not significant. All these observations indicate 
that quasiparticle decay through transitions into the bulk 
continuum is likely to play a minor role in the anisotropy 
of the scattering rates. 

It is often stated that an ideal Dirac cone prevents 
electrons from backscattering. However, as the constant- 
energy contours inside the bulk band gap become more 
and more circular with increasing BE, in Figs. if) and 
l^c) we observe a pronounced increase of the scattering 
rates along EM as compared to EK direction. This ef¬ 
fect is due to the fact that increasing the BE inside the 
bulk band gap leads to a higher probability for small- 
angle scattering events because the variation of the spin 
direction in adjacent k-points is smaller in a circle or a 
hexagon than at the tips of a snowflake-like constant- 
energy contour, as discussed above. The overall effect 


is consistent with our explanation that spin-dependent 
scattering causes the observed anisotropy of the scatter¬ 
ing rates in Bi 2 Te 3 . 

To summarize, we have analyzed the effect of strong 
hexagonal warping on the lifetime broadening of the TSS 
of the prototype TI Bi 2 Te 3 by means of angle-resolved 
photoemission. We have shown that the large Fermi- 
surface distortion introduced by warping leads to an 
anisotropy of the scattering rates which is larger along 
the EK high-symmetry direction. Based on the theoret¬ 
ically predicted behavior of the ground-state spin tex¬ 
ture of a strongly warped Dirac cone, we have identified 
spin-dependent scattering as the underlying mechanism 
giving rise to anisotropic scattering rates of the surface- 
state electrons. These results could help paving the way 
for controlling surface-scattering processes via external 
gate voltages in future spintronic devices based on TIs. 
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